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A SEMIEMPIRICAL PROCEDURE FOR ESTIMATING WING BUFFET 
LOADS IN THE TRANSONIC REGION 
By T. H. Skop inski and Wilber B. Huston 

SIWMARY 


Wing buffeting data obtained from a flight loads investigation of 
two jet-propelled swept-wing airplanes were analyzed together with some 
published wind-t\mnel results on two \mswept-wing models to study the 
effects of Mach number and reduced frequency on the wing buffet lohds. 

As a result of this analysis, a semiempirical procedure for estimating 
wing biiffet loads in the transonic speed region was derived which appears 
to correlate the data from wings of different thickness ratios, stiffness, 
and size in the Mach number range from O.65 to 1.0. 

A study of a buffet -intensity-parameter variation with reduced fre- 
quency Indicated consistent trends between the flight results and the 
wind-t'unnel results. Including penetration and wing thickness ratio as 
parameters in an expression for the root-mean-square buffet moment at 
the wing root reveals Mach number effect as being a systematic departure 
from a faired curve which represented the wind-tunnel buffet-intensity- 
parameter variation with reduced frequency. This departure appeared to 
be fixed by the transition to critical flow and was a function of the 
product of the fundamental frequency and mean aerodynamic chord of the 
wing. 


INTRODUCTION 


The concept of buffeting as the response of a simple elastic system 
to random forces associated with separated flow has been explored in a 
mjmber of analytical and experimental studies (refs . 1 to 6) . These 
analyses have indicated the validity of applying the methods of gen- 
eralized harmonic anadysis to the study of wing-buffet loads. For 
example. Information presented in reference 5 has Indicated a wind-tunnel 
technique for scaling the buffet loads obtained on simple wind-tunnel 
models to the loads expected for the airplane. 
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The purpose of the present investigation was to study the flight 
buffet-loads data obtained on two sw 6 pt-wing jet-propelled airplanes and 
the wind-tunnel restilts presented in reference 5 to examine the effects 
of Mach number and reduced frequency on the wing-buffet loads. As a res\ilt 
of this analysis a semiempirical procedtire for a rough estimation of wing 
buffet loads which appears to correlate the data from wings of various 
thickness ratio, stiffness, and size has been derived. 


SYMBOLS 


b 

c 

c 

c(y) 

t/c 

y 

yg 

q. 

t 

At 

% 





M 


wing span, ft 
wing chord, ft 

mean aerodynamic wing chord, ft 
wing chord distribution, ft 
thickness ratio 

spanwise coordinate from wing center line 

spanwise coordinate from wing center line to strain-gage 
location 

dynamic pressure, Ib/ft^ 
time, sec 

increment in time, sec 
airplane normal-force coefficient 

incremental normal-force coefficient or penetration above 
buffet boimdary 

buffet input spectrum at a given value of reduced frequency 

effective slope of lift curve for vibrating wing under con- 
ditions of separated flow 

Mach n\miber 


V 


airspeed, ft/sec 
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BM 

wi(y) 

m(y) 

s 


32 

% 

M]_ 


incremental buffet moment, ft-lb 

deflection of wing elastic axis in first wing bending, given 
in terms of a unit tip deflection 

wing mass distribution, sl\jgs/ft 

wing area, ft^ 

b/2 


effective area in bending, 2 / c(y)w 2 ^(y) dy, sq ft 


ph/2 

effective area in bending, 2 / c(y)wQ^^(y) dy, sq ft 

pb/2 


mass of wing, 2 / ni(y) dy, sl\jgs 

'-'0 


effective mass of wing in bending, 2 
slugs ^ 0 


^b/2 


m(y)w 3 _ 2 (y) dy. 


effective moment of mass outboard of point y^, 
Pb/2 

/ (y - yg)m(y)w 2 _(y) dy, slug-feet 

yg 


O’M 

(Ol 


dimensionless structiural" factor. 


Mmi 


Mill 

f eSgSMy 


physical factor 


^ “1 |\ 


ft^-ib^/^ 


root -me an- square value of incremental buffet bending moment, 
ft-lb 

circular frequency for first symmetrical wing bending, 
radians/sec 
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V 



reduced frequency 

buffet-intensity parameter per unit penetration, 

d 


AIRPLANE AND INSTRIMENTATION 




The flight test data of the present study were obtained on the North 
American F-86A fighter airplane and on the Douglas D-558-II research air- 
plane. Flight tests on the F-86A airplane were conducted at the Ames 
Aeronautical Laboratory and those on the D-^^ 8-11 airplane, at the NACA 
High-Speed Flight Station at Edwards, Calif. Both airplanes have swept 
wings and are jet-powered. Three-view drawings of the airplanes are 
shown in figure 1 and the airplane characteristics pertinent to the 
present study are summarized in table I. 

Wing shear and bending-moment measurements were obtained from elec- 
trical wire resistance strain gages mounted on the spars near the wing 
root. For both test airplanes, standard strain-gage- calibration pro- 
cedixres (ref. 7) showed that the bridge mounted near the flange of the 
rear spar responded primarily to bending moment and was a reliable meas'ure 
of the variation of the wing -root incremental bending moment. 

The strain-gage responses were recorded on an l8-channel oscillograph 
which used galvanometer elements with a natural frequency of 100 cycles 
per second and which were damped to approximately 70 percent of critical. 
Standard NACA recording instruments synchronized by a common timer were 
used to measure all the pertinent quantities. The true Mach number for 
both airplanes was determined from a nose -boom airspeed system that meas- 
ured both static and total pressure and was calibrated over the test Mach 
mmiber range. 


DATA AVAILABLE FOR ANALYSIS 


Information on the correlation between buffeting intensity and 
penetration at various altitudes and airspeed was obtained from five 
selected flights that were a part of the two flight-research programs 
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conducted with the North American F-86 a and the Douglas D-558-H air- 
plELTies. The flights analyzed were ones in which gradual pull-ups were 
made into the buffeting region. The Mach numbers for all the runs were 
held as constant as possible until a maximiim acceleration was reached 
after which the Mach number was allowed to decrease gradually at these 
higher accelerations. 

The gradual pull-up data analyzed for the F- 86 a airplane included 
nine runs of two flights at 35^000 feet at Mach numbers from about O .60 
to about 0-95 six runs of a flight at an altitude of ^ 5^000 feet 
covering a Mach number range from about 0.50 to about 0 . 90 . The two 
flights of the D-558-H airplane which were analyzed Included one run at 
55^000 feet which covered a Mach number range from about 0.80 to about 
0.95 a't an altitude of 25^000 feet which included buffeting 

data for Mach numbers between O .65 and 0-80. 

Selected for comparison are data obtained in the Langley 2- by if- foot 
flutter research tunnel for two 4-percent-thick unswept semispan models 
with an aspect ratio of 4.0 and a taper ratio of 0.20- An angle of attack 
was selected that was considered to be representative of a penetration of 
approximately O. 5 O in lift coefficient for a buffet boiindary. that was 
fairly flat. The tunnel tests were limited to Mach numbers below 
about 0-75- 

In order to indicate the general character of the flight records, 
some representative time histories of Mach number, airplane normal-force 
coefficient, incremental buffeting moment, and the root-mean- square buf- 
feting moment are shown in figure 2. The Mach number and normal- force 
coefficient were calculated every 0.10 second whereas the strain-gage 
records were read every 0.01 second in order to siarvey the loads from 0 
to 50 cycle per second. 


In order to obtain the root mean square of the Incremental buffet 
moment, the component of the wing bending moment due to the airplane 
maneuver had to be removed. This was accomplished by filtering out the 
low-frequency component due to the maneuver with a numerical filter, the 
details of which are presented in reference 6. After a time history of 
the incremental birffet moment BM(t) was obtained, it was then possible 
to break this component into a series of overlapping l/2-second intervals, 
to determine the root mean square for each interval by the fomula 




T t+.24 

50 ' 


1/2 


nnH to use It as a measure of Intensity. 
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RESULTS MD DISCUSSION 


Analysis of the bxiffet data of the present investigation has been 
guided by the analytical ej^resslon 



( 1 ) 


presented in reference 5 "the root-mean- square buffet moment at the 
\d.ng root. The first term in the brackets includes the principal physical 
characteristics of the wing. The operating conditions are represented by 
the square root of the dynamic pressure. The quantity Ps is a dimension- 
less structural factor which reduces the actual values of wing area, mass, 
and moment , arm to their effective values in the vibrating system. (See 
section on "Symbols" for expression defining these effective quantities.) 
The second bracketed term in equation (l) represents the ratio of a buf- 
fet input spectrum to the coefficient of aerodynamic damping. This ratio 
is for a particiilar value of reduced frequency corresponding to the fun- 
damental bending frequency, the average chord, and the flight speed. Thus, 
equation (l) indicates that speed changes wo\ild affect btiffet loads in 
two ways. A direct effect is through variations in dynamic pressure, the 
loads being directly proportional to the square root of the dynamic pres- 
sure and thus to the first power of the speed for a system with aero- 
dynamic damping. The loads would also be affected through changes in 
the ratio of the input forces and the damping which were functions of 
reduced frequency. 

In addition to the two speed effects noted in equation (l), buffet 
loaxis at transonic speeds have been found in references ij-, 5 ^ 9 -nd 6 to 
vary linearly with penetration above the buffet boundary so that the 
expression for the root-mean-square buffet moment at the wing root for 
transonic speed can be taken as 
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or for simplicity 


= ks^Fsd) 



(5) 


where the principal physical characteristics of the wing are represented 
by the factor 


kg = 0^ IJcSMy 


and where the aerodynamic forces are represented by a buffet- intensity 
parameter 



Wind-tunnel tests which were run on a family of related 4 -per cent -thick 
unswept wings have shown* a variation of the buffeting parameter 



with reduced frequency and have also indicated (ref. 5) ^ effect of Mach 
number at the lower values of reduced frequency corresponding to Mach 
numbers greater than about O.65. In the present analysis of the avail- 
able buffet data on two swept wings, an attempt has been made to separate 
the effects of penetration and dynamic pressure in order to examine the 
effects of Mach nimiber and reduced frequency on the buffet-intensity 
parameter. 


Variation With Penetration and Dynamic Pressure 

In order to analyze the variation of buffet intensity with penetra- 
tion for the Mach number range covered by the tests, the buffet boundaries 
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of the test airplanes vere determined. The buffet boundaries for both 
airplanes are shown in figure 3 * From a time history of M and 

the value of penetration above the buffet boundary was obtained 

for the same time that (%i(t) was calcxilated. The resulting plot of 
oj^/^ as a f -unction of penetration for different Mach numbers is shown 

in figure k. Variation of buffet intensity with altitude appears to be 
accounted for when is di-vided by ^ in agreement with analytical 

results for the case where aerodynamic damping predominates. The results 
for both airplanes indicate a reasonably linear variation of buffet inten- 
sity with penetration beyond ‘the buffet boundary for the transonic speed 
range covered, and straight lines have been faired from the origin through 
the data for each Mach number. 


Variation With Mach Number and Reduced Frequency 
In order to study the effects of Mach number and reduced frequency 

on the buffet intensity parameter values of the slope corre- 

sponding to the straight-line fairings of fig-ure V have been used together 
with values of k3 and Fg for the two wings. The parameters used to 
calculate the constant ks and the dimensionless structural factor Fg, 

are given in table I for each of the two -wings. The assumed shape of 
the deflection of the wing elastic axis in first wing bending used to 
calculate the effective areas, masses, and moments in the present study 
is given by the expression 


wi(y) = 1 - cos jty/b - 2 ^ y 2 


where the deflection is given in terms of a unit tip deflection. Cal- 
culations made by using slightly different first bending mode shapes 
(obtained from refs. 8 and 9) indicate no appreciable effect of mode 
shape on the structural factor Fg. 

The biiffeting-intensity parameter calculated by the following 
expression: 


<5 



O’M 
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for the two test alrpleines is shown in figure 5 for different values of 
reduced frequency a^c/'V. Also shown for comparison in figure 5 "the 

wind-tunnel results given in reference 5 for two unswept A-percent-thick 
wing models. The mean line faired through the model data appeaxs to 
represent a characteristic variation of the buffet-intensity parameter 
with reduced frequency. Dotted lines faired through both the flight and 
wind-timnel data appear to represent a characteristic Mach number effect, 
a reduction of intensity with increasing speed. This variation attributed 
to Mach number effect can be thoi;ight of as being associated with a change 
in the area exposed to the fluctuating pressures in separated flow, thus, 
as the shock wave moves toward the wing trailing edge with increasing 
Mach number, the exposed area, and likewise the input, become smaller. 

The extrapolation shown in figure 5 to an intercept at a reduced frequency 
corresponding to sonic velocity was made on this basis. 

Examination of the data in figure 5 shows much h igher value s of the 
intensity parameter for the thicker wing s, a result which is not in con- 
flict with the results of two-dimensional airfoil studies such as ref- 
erence 10 . This observation led to the empirical attempt to include the 
effects of airfoil thickness on buffeting intensity by simply dividing 
the intensity-parameter values presented in figure 5 by the thickness 
ratio or 


o’M = 



( 5 ) 


The resulting thickness -modified buffeting-intensity-parameter values 
are shown in figure 6 . Since the choice for the two swept wings was 
arbitrary, the root -thickness ratio perpendicular to the sweep axis was 
selected. The results in figure 6 indicate both agreement with and sys- 
tematic departure from the solid line representing the intensity varia- 
tion with reduced frequency. The reduced frequency at which this depar- 
ture occurs is the value corresponding to a Mach number of approximately 
0.65 indicated by the solid symbols. The departure from this line or 
shape is perhaps fixed by the transition to critical flow and is thus 
a function of the product coqc. 

In view of the marked differences in plan form and thickness ratio 
of the wings and the models of the present study, the consistency of the 
two sets of data shown in figure 6 is s\arprlslng and perhaps fortuitous. 
In particular, the indicated unimportance' under separated flow conditions 
of sweep, taper, and aspect ratio as compared with thickness is in marked 
contrast to the state of affairs for potential flow and verification by 
a sirLtable series of systematic tests Including higher values of sweep 
angle would be desirable. Yet the consistency shown in figiare 6 is such 
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as to encoijrage a hope that these data together with equation ( 5 ) would 
permit a roiagh estimate in the design stages of the btiffet loads for a 
given wing in the Mach number range O .65 to 1.0. For such application, 
it is believed that the linear dependence on thickness ratio should not 
be expected to hold for thickness ratios less than il- percent. 


Estimation of Wing Buffet Loads 


In order to illustrate the empirical procedure outlined in this 
paper the variation of buffet load with Mach number has been estimated 
for the If-percent- and the 10-percent-thick wings of the X-1 series air- 
planes. The areas of these two wings are the same but, in addition to 
the difference in wing thickness ratio, the wing natural frequency, mean 
aerodynamic chord, aspect ratio, and mass differ; thus, a change of all 
the parameters in equation ( 5 ) is involved. Since a limited amount of 
experimental wing-bviffet-load data (ref. 11) were available for the 
10-percent thick wing at an altitude of approximately 20,000 feet, esti- 
mates of the variation per \inlt tliickness of the buffet intensity param- 
eter with reduced frequency were made for this altitude. Figure 7 is 
presented to show the estimated reduced-frequency variation of the inten- 
sity parameter which departs from the basic curve at a value of reduced 
frequency corresponding to a Mach number of O .65 and then decreases rapidly 
to zero at a reduced frequency corresponding to sonic velocity. After 


n— J 

the estimated values of — together with the necessary parameters 

(t/c) 

in equation ( 5 ) were known, it was possible to calc\ilate the variation 
of root-mean-square buffet moment per unit penetration with Mach number. 
The calculated values of om/^^sCji shown in figinre 8 for the two wings 


indicate a much smaller wing root buffet moment for the thinner wing at 
a given penetration into the bviffet region. On the basis of the effective 
areas, masses, and moments calcoilated for the two wings, the physical 
factor k 0 influences the magnitude of the buffet load more than does 

the structural factor Fs. 


In order to obtain some indication of the accuracy of the estimated 
buffet loads for the 10-percent-thick X-1 wing the published flight resvilts 
of the X-1 airplane were analyzed. The wing-root buffeting shear loads 
presented in reference 11 are one-half the maximijm peak-to-peak shear 
fluctuations in time intervals of approximately 1/2 second. From ref- 
erence Ij- the root mean square for a l/2-second sample and a wing natural 
frequency of 14 cycles per second is estimated to be approximately one- 
half the maximum value. An assumption that the wi ng lateral centers of 
lift presented in reference 11 were approximately e qual to. the centers 
of the bi af feting load was made to obtain a value of the root-mean-square 
buffet moment. The data for Mach numbers of approximately O.'jC and 0.79 
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at an airplane lift coefficient of approximately 0.82 give the values of 
shown as symbols in figure 8. The flight-test results agree 

reasonably well with the semiempirical proced\ire for obtaining an estimate 
of the wing buffet moment. 


CONCLUDING REMARKS 


The semiempirical procedure for estimating wing buffet loads sug- 
gested by the present analysis appears to correlate the data from wings 
of various thickness ratio, stiffness, and size. It is Intended to pro- 
vide a simple engineering method of obtaining a rough estimate in the 
design stages of the "buffet loads for a given wijng in the Mach number 
range from O.65 to 1.0. 

The following remarks may be made with regards to the analysis of 
the available wing buffet data to study the effects of Mach niimber and 
reduced frequency on the wing buffet loads. Including penetration and 
wing thickness ratio as parameters in an expression for the root-mean- 
square bixffet moment at the wing root revealed Mach nijmber effect as 
being a systematic departure from a faired curve which appeared to repre- 
sent a characteristic variation of the buffet intensity parameter with 
reduced frequency. The depart\jre from the c\irve representing reduced- 
frequency effect was perhaps fixed by the transition to critical flow 
and was a function of the product of the first bending frequency and 
average chord o:^c. In view of the marked differences in plan form and 

thickness ratio of the wings and the models of the present study, the 
consistency of the two sets of data was surprising and perhaps fortuitous. 
Verification by a suitable series of systematic tests including higher 
values of sweep angle would be desirable. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 19, 195^. 
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TABLE I.- PHySICAL CHARACTERISTICS OF TEST AIRPLAMES 



Airplane 

Item 

North American 

P-86a 

Douglas 

D-558-H 

Root airfoil section 

NACA 0012-64 (normal to 

NACA 63-010 (normal to 


0.25- chord line) 

O.3O- chord line) 

Tip airfoil section 

NACA 0011-64 (normal to 

NACA 63-012 (normal to 


0.25- chord line) 

0 . 30- chord line) 

Dihedral angle, deg 

3.00 

-3.00 

Geometric twist, deg 

2.00 

0 

Sweep, deg 

5 lf .25 

55-00 


(0. 25-chord line) 

(0.30 chord line) 

Taper ratio 

0.51 

0.56 

Aspect ratio 

^.79 

5.57 

Span b, ft 

57.12 

25.00 

Mean aerodynamic chord, c f t . . . • 

8.08 

7.27 

Total area, S, sq ft 

Primary bending frequency 

287.90 

175.00 

cps 

Effective area in bending, Sx, 

8.00 

12.50 

sq ft 

Effective area in bending, S2, 

88.20 

55.00 

sq ft 

51.80 

55.00 

Wing mass. My, slugs 

w.iif 

35.40 

Effective mass, Mj, slugs 

Effective moment of mass, 

8.23 

6.52 

slug-ft 

Dimensionless structural 

146.27 

69.10 

factor, Ps 

0.18 

0.17 

Physical factor, kg, ft^-lb^/^ . . . 

31.5 X 10 ^ 

20.8 X 10 ^ 
















ra™ 


Symbol 

Airplaxie 

Sweep, 

deg 

“1C, 

ft/sec 

O 

Model 

0 

318 

O 

Model 

0 

424 

A 

F-86A 

34.2 

406 

o 

D-558-n 

35.0 

571 


icr-cr 


Mean line 
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□ Wing model 
O Wing model 
A F-86A airplane 
O D-558-n airplane 


Solid symbols for M = 0.65 














































